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We present results of numerical calculations showing a three-body orbit’s period’s T dependence 
on its topology. This dependence is a simple linear one, when expressed in terms of appropriate 
variables, suggesting an exact mathematical law. This is the first known relation between topological 
and kinematical properties of three-body systems. We have used these results to predict the periods 
of several sets of as yet undiscovered orbits, but the relation also indicates that the number of 
periodic three-body orbits with periods shorter than any finite number is countable. 
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I. INTRODUCTION 

There is, at present, no deeper understanding of pe¬ 
riodic three-body orbits in Newtonian gravity, than the 
simple change of scale of spatial and temporal coordi¬ 
nates, see Sect. 10 of Ref. [U , that can be compared with 
Kepler’s third law for two-body motion, Ref. [1]. Kepler 
extracted his laws from the astronomical data concerning 
two-body periodic orbits collected by Tycho Brahe and 
his predecessors. 

Unlike a two-body orbit, a periodic three-body or¬ 
bit is characterized both by its kinematic and geometric 
properties and by its topology, which can be described 
algebraically by a word, or an element w(a, b, A, B) of 
free group F 2 on two letters a, b (and their inverses 
A = a _1 ,B = b” 1 ), see Refs. @-0|- The algorithm used 
for “reading” of words correspondingto periodic orbits 
is described in the Appendix of Ref. [4(. 

Graphically, this amounts to classifying closed curves 
according to their “topologies” in a plane with two punc¬ 
tures. The closed curves are stereographic projections of 
periodic orbits from the shape-sphere, with three punc¬ 
tures - for a detailed explanation, see Refs. M , and for 
graphic illustrations, see the web-site [lf| - onto a plane 
with two punctures, the puncture at the “north pole” 
having been projected to infinity. 

That procedure leads to the aforementioned free group 
F 2 on two letters (a, b), where (for definiteness) a de¬ 
notes a clockwise “full turn” around the right-hand-side 
puncture, and b denotes the counter-clockwise full turn 
around the other puncture in the plane/sphere. In this 
way the topology of an orbit can be transformed into an 
algebraic object that can be further manipulated. 

But, even within this particular method of assigning a 
sequence of symbols to a topology there remains an am¬ 
biguity, regarding the question which puncture should 
be taken as the “north pole” of the stereographic projec¬ 
tion, see Appendix I A II The length of the word generally 
depends on this choice, see Appendix IB II We resolve 
this ambiguity by using the (common) symmetry axis of 
all presently known collisionless zero-angular-momentum 


periodic orbits to define the “north pole”, which we call 
the “natural”, or “symmetric” choice, because it leads to 
equal numbers n a = rig = \n w of small letters a and b, as 
well as equal numbers tia = riB = \n w of capital letters 
A, or B. These relations need not hold with a different 
choice of “north pole”, however, e.g. with cyclically per¬ 
muted punctures, generally m a ^ nib and niA ^ nie, see 
Appendix IB II Moreover, the above-described procedure, 
is not the only way of assigning a sequence of symbols to 
a topology, for an alternative, together with our results 
expressed in these alternative terms, see Appendix IB 21 

To date there is no collection of astronomical data re¬ 
garding periodic orbits of three bodies comparable to 
Brahe’s collection of two-body orbits. Therefore, if one 
wishes to study general properties of the three-body sys¬ 
tem one must resort to numerical studies. To this end, in 
this Letter we use the world’s total (published) data set 
containing 46 distinct collisionless periodic orbits, Refs. 
mm , to extract the following (wholly unexpected) 
linear dependence of the (generalized) Kepler’s third law 
for the ratio T r (w)/T r (w$) of “rescaled” periods (i.e. pe¬ 
riods evaluated at one common energy E ) of three-body 
orbits, 

T r (w) _ T e (w)\E{w)\ 3/2 ^ N w _(n w +n w ) 

T r (w s) ~ T e (w s )\E(w 8 )\ 3 / 2 2 2 U 

on their topologies w, specifically on (one half of) the 
number of all letters N w = (n w + n w ), see Fig. |T] Here 
n w is the number of small letters a, or b, and n w is the 
number of capital letters A, or B contained in the latter 
w, and w 8 = abAB if the free-group word describing the 
figure-8 orbit, Ref. ill- We have divided the total 
number of letters N w into two parts because orbits fall 
into different classes with distinct values of n w and n w , 
see Table HI 


The worst-case disagreement of this linear dependence 
with the data is about 10%, though generally it is (much) 
better, occasionally reaching the limits (six significant 
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Figure 1: (color on line) The rescaled periods T r (w) of 
presently known (zero-angular-momentum) three-body orbits 
divided by the period of the figure-8 orbit T r (w s), versus one 
half of the length of word N w , i.e., one half of the number 
of all letters in the free-group word w describing the orbit, 
N w /2 = (n w + n w )/2, where n w is the number of small letters 
a, or b, and n w is the number of capital letters A, or B in 
the letter w. Four (linear, quadratic, cubic and quartic) fits 
are shown as solid lines of different colors, yet they overlap so 
much that the difference can be seen only at N w > 80. 

decimal places) of our numerical precision. In this Letter 
we point out four more accurate (than 10 %), yet still 
approximate sub-sequences of orbits, and one possibly 
exact regularity. For clarity’s sake, we show in Fig. [2] the 
corresponding graph for non-satellite orbits, only. Note 
that the range of the abscissa in Fig. [2] only reaches the 
value N w = 49. 


At this point a few words must be said about the statis¬ 
tical significance of results presented in Fig. [I] or equiv¬ 
alently about the number of distinct periodic 

orbits of “length” N w : 1) At small length values one 
can explicitly count the mathematically allowed orbits 
and show that many have already been found, see Ap¬ 
pendix [A~2l 2) As one increases the number of letters N w , 
the number of topologically distinct orbits #(N W ) grows 
rapidly, see Appendix IA 21 and the number of presently 
discovered (and displayed here) orbits pales in compari¬ 
son with that number. The number jfc(N w ) is not nec¬ 
essarily the same one as the number of physically possi¬ 
ble orbits - Moore has shown by explicit examples how 
mathematically allowed orbits disappear as the exponent 
in the potential is reduced from a = 2 to a = 1 in New¬ 
tonian gravity, Ref. 0 

The large number of still possibly undetected orbits 
makes the observed linearity of the graph, Fig. [Q at 
higher values of N w all the more impressive: Note that 
24, out of grand total of 46 orbits taken from Refs. 
0, H 0 0, extend up to N w = 49. These 24 orbits in¬ 
clude 10 (non-choreographic) figure-eight satellites from 
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Figure 2: (color on line) The rescaled periods T r (w) of 16 
presently known non-satellite zero-angular-momentum three- 
body orbits divided by the period of the figure-8 orbit T r (ws) 
versus one half of the length of word N w , i.e., one half of the 
number of all letters in the free-group word w describing the 
orbit, N w /2 = (n w + n w )/2, where n w is the number of small 
letters a, or b, and n w is the number of capital letters A, or 
B in the letter w. Four (linear, quadratic, cubic and quartic) 
fits are shown as solid lines of different colors. 


Ref. dj. Among these 24 there are 16 non-satellite orbits 
that are shown separately in Fig. [2] The remaining 22 
(of 46) orbits are the new (fc=5, 7, 14, 17, 22, 26, 35, 
41 figure-eight satellite) choreographies, Ref. that 
extend up to N w = 82 and thus test the proposed lin¬ 
ear dependence(s) farther into the previously unexplored 
region. We emphasize that three of these new choreo¬ 
graphic orbits are not satellites of the figure-eight. 

The most precise regularity explains two previously no¬ 
ticed, but unexplained identities: a) the identity of peri¬ 
ods (to 16 decimal places, in Ref. Q); and b) the identity 
of actions (to seven significant digits, in Ref. 0) , of two 
distinct orbits with the same topology, viz. of Moore’s 
and Simo’s figure-8 solutions, evaluated at equal energies. 
The same phenomenon was observed among seven inde¬ 
pendent k = 17 figure-eight satellites reported in Ref. 
0, and the two butterfly orbits in Ref. [ll|. Moreover, 
the ratio of periods of two different-fc satellites equals the 
ratio of the corresponding k' s, to high precision. 

An extension of this rule to arbitrary k predicts the 
periods, and actions of new, as yet undiscovered three- 
body orbits and helps one to identify and classify newly 
discovered orbits, as in Ref. @. The four less precise 
linear regularities also predict new families of orbits, with 
somewhat lower precision. Before we proceed to present 
the evidence and state the regularities in their algebraic 
forms, we must first define what we mean by Kepler’s 
third law in a three-body system? 
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Kepler’s laws for three bodies? 

Kepler’s first law, which states that a periodic two- 
body trajectory must be an ellipse, manifestly holds only 
in the two-body case; his second law, which asserts con¬ 
servation of angular momentum, holds for any number 
of bodies, but generally does not have a simple geomet¬ 
rical interpretation as in the two-body case. His third 
law, which relates the square of the orbital period T of a 
planet to the cube of the semi-major axis a of its orbit, 
i.e., T 2 oc a 3 , also has a proper three-body generalization, 
despite its appearance. As only Lagrange’s three-body 
orbits are elliptical, we must address the question what 
corresponds to the semi-major axis a in Kepler’s third 
law for arbitrary periodic three-body orbits? 

The answer to this question can be found by remem¬ 
bering the virial theorem, Ref. [l|, for the Newtonian po¬ 
tential and then by using the hyper-spherical three-body 
variables, Ref. Q- As the “overall size” of a three-body 
system is determined by a time-average (to be specified 
below) over one period of the hyper-radius R ; and the 
virial theorem relates the time-average of the potential 
energy V(R ) ~ 1/7? to the (constant) total energy E, so 
it is the harmonic mean (over one period) of the hyper¬ 
radius R , defined as jj = y f 0 that corresponds to 
the semi-major axis a in the three-body case, see footnote 
& 

Thus, we (may) replace the “mean size” R of the three- 
body system in Kepler’s third law T oc R 3 / 2 with the 
inverse absolute value of energy 17?| 3 , i.e., T oc |i?| _3//2 , 
or equivalently T|.E| 3 / 2 = const. [22|. The constant on 
the right-hand-side of this equation is not “universal” in 
the three-body case, as it is in the two-body case: it may 
depend on both the family of the three-body orbit and 
its angular momentum L, see Refs. [3, [3 • The angular 
momentum L transforms differently under scale transfor¬ 
mations than the period T, which is the reason why only 
the vanishing angular momentum L = 0 is a “fixed point” 
under scaling. For this reason, we confine ourselves here 
to three-body orbits with zero angular momentum, and, 
for simplicity’s sake, to equal masses. 


II. OBSERVED KEPLER-LIKE REGULARITIES 

As stated in the Introduction, families of periodic 
three-body orbits can be characterized by the topol¬ 
ogy of their trajectories in the real configuration space 
(“braid group”), or on the (so-called) shape sphere (“free 
group”), as described in Refs. [2H5j|, the latter can be 
specified by the conjugacy classes of elements, or “words” 
w , for short, consisting of letters a, b, A = a'/B = b _1 , 
that define the free group on two letters (a, b). 

Thus, here we must study the dependence of the 
constant on the right-hand-side of the scaling law 
T(w)\E(w )\ 3//2 = const (w) on the structure of the word 
w(a, b, A, B) that characterizes a periodic three-body or¬ 


bit with zero-angular-momentum. 

Numerical observations 

Note that if we classify orbits using only the num¬ 
bers, n = n a = rib and h = n* = of a,b and A,B 
letters respectively, in the orbits’ free group group ele¬ 
ments, then all periodic orbits satisfy roughly the linear 
rule T ~ n + fi, c.f. column and (n + n)/ 2, in Table 
ID and Fig. [Q By perusal of Table Q] one finds another 

Table I: Rescaled periods T of three-body orbits, their ra¬ 
tios with Moore’s figure-8 period Tms, and with period 
of the first orbit /3 in the given section of the Table, as func¬ 
tions of the numbers n^,n b ,nA,n B , of a’s, b’s, A = a -1 ’s and 
B = b _1 ’s respectively, in the free-group word description of 
the orbit. Note that “generic” relations n = n a = n b and 
n = 71a = riB hold only in the “natural” or symmetric choice of 
stereographic projection (see the Introduction and Appendix 
lAlT ) and that, due to the time-reversal symmetry of the solu¬ 
tions, the n and h may be interchanged. 


Label 

<T>2 

<T>2 

<t> 2 

n-\-n 

(n,n) 

(Tm8)2 

<T B >2 

n0+fi0 

M8 

26.1281 

1 

1 

l 

1,1 

S8 

26.1268 0.999951 

0.999951 

l 

1,1 

I.B.l moth I 

68.4636 

2.62031 

1 

l 

2,3 

II.B.l yarn 

205.469 

7.86391 

3.00114 

3 

6,9 

I.A.l butterfly I 

56.3776 

2.15774 

1 

1 

2,2 

I.A.2 butterfly II 

56.3746 

2.15762 

0.999944 

1 

2,2 

I.B.5 goggles 

112.129 

4.29152 

1.9889 

2 

4,4 

I.B.7 dragonfly 

104.005 

3.98059 

1 

1 

4,4 

I.A.3 bumblebee 

286.192 

10.9534 

2.7517 

11/4 

11,11 

II.C.2a yin-yang I 

83.7273 

3.20449 

1 

1 

3,3 

II.C.2b yin-yang I 

83.7273 

3.20449 

1 

1 

3,3 

II.C.3a yin-yang II 334.876 

12.8167 

3.9996 

4 

12,12 

II.C.3b yin-yang II 334.873 

12.8166 

3.9996 

4 

12,12 

I.B.l moth I 

68.4636 

2.62031 

1 

1 

2,3 

I.B.3 butterfly III 

98.4354 

3.76742 

1.43778 

7/5 

3,4 

I.B.2 moth II 

121.006 

4.63126 

1.76745 

9/5 

4,5 

I.B.4 moth III 

152.33 

5.83013 

2.22498 

11/5 

5,6 

I.B.6 butterfly IV 

690.627 

26.4324 

10.0875 

49/5 

24,25 


five sequences of orbits, separated by horizontal lines in 
Table [T] The first one of the five sequences in Table Q] 
contains two orbits that form in a progenitor-satellite 
relationship: w(II.B.l yarn) = w(I.B.l MothI) 3 , where 
the “II.B.l yarn” solution whose topology is described 
by w(II.B.l yarn) = (ba(BAB)ab(ABA)) 3 , i.e., as the third 
power (jfe = 3) of the (time-reversed) “I.B.l Moth I” orbit 
described by u>(I.B.l MothI) = ba(BAB)ab(ABA). We can 
see in Table U that the ratio of their periods equals three, 
to one part in 3 x 10 4 . This constitutes the presently 
available data about satellites of orbits other than the 
figure-eight one. The remaining four sequences in Table 
U define “empirical integer relations/laws” with linearly 
related periods - we record them here, both as a challenge 
to mathematicians, and as a guide to future numerical 
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searches, as they predict periods of as yet undiscovered 
orbits. 


Predictions 

On the basis of these empirical regularities, we predict: 

1. new “yin-yang” orbits with ratios of periods 

T (yyk)/ r (yyi) = k = 2 , 3 , 5 ,.. 

2. new “butterfly I - goggles” orbits with ratios of pe¬ 
riods T/T(I.A.l) equal to 3,4,5, ... ; 

3. new “dragonfly - bumblebee” orbits with 
T/T(I.A.3) = 5/4,3/2, 7/4,.. 

4. new “moth I,II,II - butterfly III” orbits with n = 
6, 7,.... The “butterfly IV” orbit deviates the most 
(8%) from this sequence, and may well define a sub¬ 
sequence of its own. 

5. new satellites of “moth I”, above and beyond the 
“yarn” orbit. 

III. THREE POTENTIALLY EXACT 
PROPOSITIONS 

The above numerical observations led us to search for 
and examine other examples of integer powers of orbits, 
such as the satellites of the figure-eight solution, Ref. 
0 , 1 - 

Statement 

Thus, we found the following regularities/laws in Table 
HI1 that hold to a higher numerical precision than T ~ 
n + n: 

1. All periodic orbits, described by one topology w 
and normalized to a common energy E, have almost 
identical (i.e. valid to a precision comparable to the 
numerical precision of the solutions) values of the 
period Ti(w) = T 2 (w) = T 3 (w) = ... and of the 
(minimized) action S\(w) = S 2 (w) = S 3 (w) = ■■■', 

2. Kepler’s third law “constant” T{w k )\E(w k )\ 3/2 of 
the “fc-th satellite orbit” (specified by the free- 
group element w k where fc is an integer) of the “pro¬ 
genitor orbit” w equals k times the Kepler’s third 
law constant T(w)\E(w)\ 3 ^ 2 of the progenitor orbit 
w: T(w k )\E(w k )\ 3 ^ 2 = kT{w)\E{w)\^' 2 . More sim¬ 
ply, periodic orbits normalized to a common energy 
E have periods related by T(w ki ) = kiT(w) [22|. 

3. all periodic orbits described by the topologies w ki , 
where ki = 1,2,3... have (minimized) actions 
S(w ki ) that are directly proportional to the expo¬ 
nent ki : S(w ki )\E(w ki )\ 1 ^ 2 = kiS(w)\E(w)\ 1 / 2 . 


Of course, proposition 1. is a special case (ki = 1) of 
propositions 2. and 3.. The fact that proposition 2. is 
identical to proposition 3. for Newtonian gravity follows 
from the identity S = 3\E\T, where S is the (minimized) 
action of a periodic motion, E is its energy and T its 
period. In the following we show that proposition 2. is 
true - by numerical examples; an outline of an argument 
regarding proposition 3., based on Cauchy’s residue theo¬ 
rem and the periodicity of three-body motion in the case 
of Newtonian gravity potential, can be found in Ref. 0] . 

Numerical evidence 

We use the numerical data from Refs. 0 , 0 : periods of 
orbits (with unit masses rrij = 1, i = 1,2,3 and Newto¬ 
nian coupling constant G = 1) from Ref. 0 and rescale 
it to the common energy E = —1/2, which rescaled pe¬ 
riods are tabulated in Table |U Moreover, we rescale the 
periods of “satellites of figure-eight” orbits from Ref. 0 
in the same way, and tabulate them in Table in 23 . 

Topologies of “satellite of the figure-eight” orbits’ are 
the fc-th powers of the figure-eight orbit (formally, their 
homotopy class is (abAB) fc ). The initial conditions for 
Moore’s “figure-eight” choreography and Simo’s figure- 
eight orbit are labelled by F 8 and 58, respectively, in 
Table EH Next note that the ratio of the period T and 
the period of Moore’s figure-8 orbit Tg in Table Ell equals 
the “slalom/satellite power” fc = T/Tg to (at least) four 
significant decimal places. Moreover, as noted by Galan 
et al. 0 , the actions of these two orbits are also identi¬ 
cal, again, up to numerical precision, in accordance with 
our proposition 1). Other examples of identity of periods 
(and actions) of two, or more distinct orbits from the 
same topological sector can be found in Table EH (see 
sectors fc = 7,11,17 and comments in Ref. 0 about 
the identity of various orbits) and in Table El see orbits 
I.A.l butterfly I and I.A.2 butterfly II. These form the 
numerical evidence for our proposition 1.. 

Table EH also presents the first numerical evidence for 
propositions 2. and 3. “topological scaling law”, albeit 
limited to the set of satellites of figure-eight solutions. 
Propositions 2. and 3. have been confirmed, after the 
fact, in a new set of figure-eight satellite choreographies 
with fc € [5,41], as described in Ref. [g], where they 
greatly facilitated the detection of three choreographic 
orbits that are not satellites of the figure-eight. 


Discussion 

1) The ratios of the shortest periods T(/3) in each of 
the four sequences of orbits in Table El and the figure- 
eight orbit Xm 8 are all larger than, or equal (“dragon¬ 
fly” ) to unity; thus it appears that the figure-8 orbit has 
the shortest scale-invariant period T\E\ 3 / 2 of all non¬ 
colliding zero-angular momentum orbits. This means 
that a search for periodic orbits with periods shorter than 
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Table II: Periods of three-body orbits rescaled to a common 
energy E = —1/2, together with the rescaling factor A — 
—2 E, where E is the energy of the orbit before scaling. T is 
the period, Tg is the period of Moore’s figure-8 orbit and k 
is the “slalom power” (i.e. (abAB) fc is homotopy class of the 
orbit). We also list Moore’s (M 8) and Simo’s (S 8) figure-eight 
orbits, for comparison. (T) represents the average value of 
three “measurements” made with different versions of Runge- 
Kutta4 and definitions of the period; (T )2 is the average of the 
“best pair”, i.e., eliminating one measurement that deviates 
the most from the arithmetic mean; for more on this, see Ref. 
0 - 


# a (t) 2 (t) (t)/(t M8 ) (r) 2 /(r Ma )2 fc 


M 8 2.57428 26.1281 26.1284 1 1 

S8 2.58387 26.1268 26.1294 1.00004 0.99995 


NCI 

3.0086 

182.873 

182.870 

6.99890 

6.99910 

7 

NC2 

2.16153 

182.873 

182.871 

6.99894 

6.99910 

7 

Ol 

3.14090 

182.854 

182.856 

6.99843 

6.99837 

7 

02 

1.99620 

182.854 

182.855 

6.99839 

6.99836 

7 

03 

3.06222 

287.333 

287.334 

10.9971 

10.9971 

11 

04 

2.97791 

287.334 

287.336 

10.9972 

10.9971 

11 

05 

2.98251 

365.760 

365.761 

13.9987 

13.9987 

14 

06 

3.12702 

444.154 

444.155 

16.9991 

16.9991 

17 

07 

3.12618 

444.154 

444.156 

16.9992 

16.9991 

17 

08 

3.10403 

444.155 

444.155 

16.9991 

16.9991 

17 

09 

3.10267 

444.155 

444.156 

16.9992 

16.9991 

17 

010 

3.06755 

444.156 

444.157 

16.9992 

16.9992 

17 

Oil 

3.06576 

444.156 

444.157 

16.9992 

16.9992 

17 

012 

3.04027 

444.157 

444.158 

16.9992 

16.9992 

17 

013 

2.92754 

444.167 

444.169 

16.9996 

16.9996 

17 

014 

2.89998 

444.169 

444.170 

16.9997 

16.9997 

17 

015 

2.85692 

444.169 

444.170 

16.9997 

16.9997 

17 


some finite prescribed value can only lead to a finite num¬ 
ber of distinct topologies. If we now assume that each 
topology has (only) a finite number of satellites for each 
value of exponent k, then we are assured that the number 
of all such orbits must be finite, and that the number of 
periodic orbits with arbitrarily long (including infinite, 
whatever that may mean) periods is infinite, but count¬ 
able/denumerable. 

2) The orbits related by the above four “approximate” 
sequences generally fall into different “algebraic” and/or 
“geometric classes” defined in Refs. [3, [HJ, which shows 
that the topology of an orbit, and its belonging to one of 
the four sequences, are not directly related to the classi¬ 
fications based on either of these two concepts. There¬ 
fore, we conclude that the newly-predicted orbits may fall 
outside of any previously known algebraic and geometric 
classes and that their initial conditions may lie outside 
of the currently explored subspace of initial conditions. 
We also note that all of the above orbits share a common 
“geometrical symmetry” on the shape sphere: they are 
all symmetric under reflections with respect to the axis 
passing through the center of the shape sphere and the 
Euler point defining the initial conditions. 

3) As noted in Refs. SI Montgomery’s jlj prescrip¬ 


tion for assigning a free-group element conjugacy class 
to a topology is not unique, see Appendix IA II Differ¬ 
ent conventions as to the North Pole of the stereographic 
projection yield different free-group words, often of differ¬ 
ent lengths. Thus, the total length of a free-group word 
describing the topology is not a constant, but generally 
depend on the convention taken for the stereographic pro¬ 
jection from the shape-sphere onto the plane with two 
punctures. 

Here we have used only the “natural”, or “symmetric” 
convention/rule, i.e. we consistently chose the puncture 
“opposite” to our initial conditions’ Euler point to be 
the North Pole of the stereographic projection. This is a 
’natural” choice because all known zero-angular momen¬ 
tum equal-mass solutions are symmetric under reflections 
about the axis defined by that Euler point and the center 
of the shape sphere, see graphs in Ref. [161] . 

This convention yields a unique free-group word for 
any given orbit with this symmetry, with a unique length. 
It is not clear, however, that all equal-mass zero-angular- 
momentum collisionless periodic orbits have to share this 
reflection symmetry. This convention can be extended to 
the case with two equal masses, because this symmetry 
“survives” one distinct mass, but not all three. 

4) After this Letter was submitted, it was pointed out 
to us that a recent mathematical preprint Ref. E2 con¬ 
tains results similar to ours, viz. orbits with periods pro¬ 
portional to the number of letters in their words, with the 
caveat that they hold only for certain newly described or¬ 
bits with non-zero angular momentum, that pass “very 
near triple collisions”. Such orbits have been obtained 
only formally, i.e. without a computed example, because 
triple collisions cannot be regularized. As these orbits 
do not satisfy two of our conditions - vanishing angular 
momentum and collisionless orbits, we disregard them in 
this context. 


IV. SUMMARY, CONCLUSIONS AND 
OUTLOOK 

In summary, we have analyzed the topological depen¬ 
dence of the scale-invariant period T r = Te\E\ 3 / 2 of 
presently known periodic three-body orbits and found re¬ 
markable (possibly exact) regularities that can be derived 
from the corresponding relation, S(w k )\E(w k )\ 1 ^ 2 = 
kS(w)\E(w )\ 1 ^ 2 , between the action S(w k ) and energy 
E(w k ) of the fc-th satellite orbit described by w k , on one 
hand, and the action S(w), and energy E(w) of the pro¬ 
genitor orbit w, on the other. 

The above proposition is, to our knowledge, the first 
published Kepler-like law for arbitrary three-body peri¬ 
odic orbits. Besides this one, possibly exact new law, 
there is empirical evidence for other approximate, but 
still unexplained regularities. 

The results presented here are striking and unexpected. 
If their validity persists after further improvements of nu¬ 
merical precision, they will be signs of deeper mathemat- 
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ical truths. But, we must at all times keep in mind the 
fact that these are the results of numerical experiments, 
which are just what their name implies: experiments. So, 
it is clear that numerical experiments cannot be used to 
prove theorems; that is left to mathematicians. 

From the practical point of view, perhaps the most 
important results are those that point towards new, as 
yet undiscovered periodic solutions. Our results do not 
predict the initial conditions of the new orbits, but only 
their periods, at given energy. 

From the conceptual point of view, perhaps the most 
important result of this study is that it points towards 
a denumerable, and perhaps even finite number of peri¬ 
odic orbits with periods shorter than some pre-specified 
finite value. If the latter (finite number) proposition is 
true, that would (greatly) constrain the number of possi¬ 
ble periodic orbits, and thus would imply that one could 
meaningfully search for and perhaps even find all periodic 
orbits with periods shorter than some prescribed value. 
At any rate, our results indicate that periodic orbits do 
not form a dense set in the set of all three-body orbits. 

These results also pose several questions in physics and 
astrophysics: 1) astrophysics: how can these three-body 
systems be formed in an astrophysical process? One 
mechanism for the formation of the figure-eight orbit 
was suggested by Heggie, Ref. [3, and it can be tried 
for other orbits. Before one attempts answering such a 
question, however, one must establish which of the new 
orbits are stable, as only stable orbits have non-vanishing 
probabilities of being formed. 2) physics: what happens 
when one replaces the Newtonian gravity with Coulom- 
bian electrostatics? In that case one ends up with new 
free parameters (ratios of electric charges), as these cou¬ 
pling constants do not dependent on the masses. Do our 
regularities survive in the case of electrostatics, and un¬ 
der which conditions? 

We conclude by noting that our results also present 
new mathematical questions, such as: 1) how many de¬ 
generate minima of the action (integral) are there in each 
topological sector? This particular question falls into the 
realm of the variational calculus in the large, or Morse 
theory, see [lj|; 2) how many algebraically distinct con- 
jugacy classes of F 2 correspond to any given topology 
of a three-body orbit? This “automorphism” question is 
an open problem in combinatorial group theory, see Ap¬ 
pendix |Xl] and Ref. J20| ■ 3) how many topologically dis¬ 
tinct orbits/algebraically distinct conjugacy classes are 
there as a function of the word length N w (see Appendix 
I A 2D ? This is a combinatorial problem. These and other 
questions remain as a challenge to mathematicians. 
Acknowledgments This work was supported by the 
Serbian Ministry of Education, Science and Technolog¬ 
ical Development under grant numbers 01 171037 and 
III 41011. We thank Dr. Mitsuru Shibayama (Osaka U.) 
for providing one set of measurements of periods and for 
a thoughtful reading of the manuscript, and Dr. Ayumu 
Sugita (Osaka City U.) for alerting us to the discrepancy 
in the measured value of figure-eight orbit’s period and 


for letting us use his code to measure other periods used 
in this paper. Last, but not least we thank Aleksandar 
Bojarov for his help with generating data in Table Hill 


Appendix A: Counting topologically distinct orbits 

1. Ambiguities in algebraic description of 
topologies 

The ambiguity mentioned in the Introduction is best 
illustrated by a simple example, already shown in Refs. 
0,1: the Broucke-Hadjidemetriou-Henon (BHH) fam¬ 
ily of orbits. This family is described by a simple loop 
around any one of three poles/punctures in the shape- 
sphere. By applying the stereographic projection of the 
shape sphere onto a plane with two punctures, a single 
loop can be “translated” into three different “words”, as 
can be seen below. 

1) If one chooses as the “north pole” of the stereo¬ 
graphic projection either one of the two punctures that 
are not encircled by the trajectory, then the orbit is de¬ 
noted by the free-group letters (a, or B). 

2) If, on the other hand, one chooses the same punc¬ 
ture that is being encircled as the “north pole”, that 
corresponds to aB, i.e., a loop around both punctures in 
the plane. 

But a single loop around any one of the three punc¬ 
tures on the original shape-sphere must be equivalent to 
a loop around either of the two remaining punctures; con¬ 
sequently, aB is equivalent to a and/or B. 

Of course, the time-reversed solutions are topologically 
equivalent to the original orbits, i.e., a = A, B = b, aB = 
bA. Thus, we see that (up to six) different conjugacy 
classes of possibly different lengths can be assigned to 
one and the same topology. 

The length of a free-group word is constant under such 
“changes of punctures” for orbits that are invariant un¬ 
der cyclic permutations of three punctures on the shape 
sphere. So far, there are only four known examples of or¬ 
bits with a constant free-group word length: one of them 
is the figure-eight orbit (and its satellites); the remaining 
three are the new non-figure-eight choreographies found 
in Ref. slj- It is an open question if any periodic three- 
body equal-mass orbits exist that have no remnant of the 
full S 3 permutation symmetry group? 

The transformations/mappings of words defined by 
rules such as those shown above are known in mathemat¬ 
ics as automorphisms of a free group - in this regard we 
can do no better than to quote Vogtmann: “The study of 
automorphism groups of free groups in itself is decidedly 
not new; these groups are basic objects in the field of 
combinatorial group theory, and have been studied since 
the very beginnings of the subject. Fundamental contri¬ 
butions were made by Jakob Nielsen starting in 1915 and 
by J.H.C. Whitehead in the 1930’s to 1950’s. ... these 
groups have attracted a tremendous amount of research 
work in spite of gaps of decades in the sequences of pa- 
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pers that deal with them and, in spite of the work that 
was done there was still much more unknown about these 
groups than known.” Ref. Jfj. In summary, the effects 
of automorphisms on conjugacy classes, and in particular 
on the number of classes equivalent under such automor¬ 
phisms, are not known in general. But, computations 
of conjugacy classes of specific words can be readily done 
using an electronic computer, see Annendix lB II One can 
also count the number of conjugacy classes of Fi explic¬ 
itly, “by hand”, for small word lengths, and/or try and 
establish upper and lower bounds on this number, as in 
Appendix I A 21 


2. Counting conjugacy classes of free group 
elements 

We are indebted to referee C for raising this question. 

1) The number of distinct AT-letter words in alphabets 
with four letters is 4A; the number of conjugacy classes 
equals 4 N /N. 

2) The above estimate of 4 N /N is a gross overestimate 
of an upper bound on the number of conjugacy classes 
in the free group F%, because one did not take into ac¬ 
count the reduction in length of all words containing two 
neighbouring identical small and large letters, e.g. (a and 
A), or (b and B). Such words immediately boil down to 
words with [N — 2) letters, etc. This reduces the number 
of conjugacy classes of F 2 (far) below A N /N\ 

3) Just how much smaller is the correct number than 
4 N /N? We cannot say precisely, but we have a lower 
bound instead: the number of conjugacy classes can not 
be smaller than (2 N /N) x const.. Taking only two small, 
or two capital letters, or one small and a different capital 
letter, such as (a and B) into consideration, the two can 
not cancel; of course there is more than one such choice, 
and that multiplicity is covered by the constant const.)! 
This constant is larger than 2, but less than 8, so it makes 
a difference only at small N. At N < 10, this number is 
generally around the square root of the upper bound. 

4) Thus we have an upper and a lower bound both of 
which grow as e dxN /N, as N — > 00 , albeit with different 
values of constant d: for the upper bound d = log 4 and 
for the lower bound d = log 2. It seems natural to assume 
that the actual number of conjugacy classes also follows 
the e de * xN /N asymptotic law, with a exact value of d ex . 
lying between the above two values: log 2 < d ex . < log 4. 

5) None of the above bounds and estimates took into 
account the automorphism ambiguity described in Ap¬ 
pendix IA II above, which means that the number of dis¬ 
tinct topologies can only be smaller than these estimates. 

6) There are further constraints, stemming from phys¬ 
ical symmetries and dynamical considerations, which 
must be imposed on the number of conjugacy classes of 
F -2 that are equivalent under automorphisms, in order for 
that number to be equal to the number of distinct orbits: 

a) Time-reversal invariance demands that an orbit and 
its time-reversed version (“inverse”) be viewed as equiv¬ 


alent. Needless to say, conjugacy classes of an (non¬ 
trivial) element and its inverse are algebraically differ¬ 
ent. Therefore, the above number of conjugacy classes 
must be halved in order to obtain the number of distinct 
orbits. 

b) Not all topologies allow zero angular momentum or¬ 
bits without collisions. The simplest example is the BHH 
family, which at zero angular momentum leads to the v. 
Schubart orbit, which is a collinear, therefore colliding, 
periodic orbit. 

The exact computation of the number of topologically 
distinct orbits as a function of word length N w is a dif¬ 
ficult combinatorial problem that remains as a challenge 
to mathematicians. 


Appendix B: Alternative symbolic sequences 
1. Alternative word readings 

We can readily implement cyclic permutations of three 
punctures in our free-group word-reading algorithm, see 
the Appendix in Ref. [J| , and this results in the new word 
lengths shown in Table IIIII Note that many of the ratios 

Table III: Rescaled ratios of periods T of three-body orbits 
and Moore’s figure -8 period Tms, and ratios with period 
of the first orbit /? in the given section of the Table, as func¬ 
tions of the numbers n a , n b , JTa, Ub, of a’s, b’s, A = a _1 ’s and 
B = b _1 ’s respectively, in the free-group word description of 
the orbit for both the “natural choice” and cyclic permu¬ 
tations of punctures on the shape sphere. Note that the 
“generic” relations n = n a = n b and n = ru = n B which hold 
with the “natural” or “symmetric” choice of punctures, do 
not hold with cyclically permuted punctures, i.e., ni a 7 ^ ni b 
and ms 7 ^ ni B . Here we define n 1 = |(ni a + m b ) and 
hi = |(n ia + ni B ). The ni and hi may still be interchanged. 


Label 

(T>2 

(T)2 

n+n 


ni+ni 


(Tm8 ) 2 

<T 3 >2 

n B +ng 

[n, ft) 

n 10~\-™10 


M 8 

1 

1 

1 

(id) 

l 

( 1 , 1 ) 

S8 

0.999951 

0.999951 

1 

(i,i) 

l 

( 1 , 1 ) 

I.B.l 

2.62031 

1 

1 

(2,3) 

l 

(2.5,3) 

II.B.l 

7.86391 

3.00114 

3 

(6,9) 

3 

(7.5,9) 

I.A.l 

2.15774 

1 

1 

( 2 , 2 ) 

1 

(2.5,2.5) 

I.A.2 

2.15762 

0.999944 

1 

( 2 , 2 ) 

1 

(2.5,2.5) 

I.B.5 

4.29152 

1.9889 

2 

(4,4) 

2 

(5,5) 

I.B.7 

3.98059 

1 

1 

(4,4) 

1 

(4,4) 

I.A.3 

10.9534 

2.7517 

11/4 

( 11 ,H) 

11/4 

( 11 , 11 ) 

II.C.2 

3.20449 

1 

1 

(3,3) 

1 

(3.5,3.5) 

II.C.3 

12.8167 

3.9996 

4 

( 12 , 12 ) 

4 

(14,14) 

I.B.l 

2.62031 

1 

1 

(2,3) 

1 

(2.5,3) 

I.B.3 

3.76742 

1.43778 

7/5 

(3,4) 

17/11 

(4,4.5) 

I.B.2 

4.63126 

1.76745 

9/5 

(4,5) 

19/11 

(4.5,5) 

I.B.4 

5.83013 

2.22498 

13/5 

(6,7) 

25/11 

(6,6.5) 

I.B .6 

26.4324 

10.0875 

49/5 

(24,25) 

123/11 

(30.5,31) 


and n "j~_ in Table IIIII are identical, but just 

nip+nip np+np ' - 1 ’ J 

as many are not: in particular no ratio in the “moth I - 



















butterfly IV” sub-sequence is independent of the choice 
of puncture. Perhaps more importantly, the ratios of 
— " 1 +" 1 —, which determine the ratios of periods , 

"lMg+niMg ’ r Jms ’ 


i.e., the slopes of linear sub-sequences are often not-as- 
good as with the “natural choice” : for example, 

for “moth I” orbit, the measured value +— = 2.62031 


is closer to —— = 2.5 than to 

riMg+iiMg 


ni+ni 

niMg+filMg 


2.75. 


Similarly, for the I.A.l “butterfly I” orbit, the measured 
value +— = 2.15774 is closer to —+21— = 2 than 

J-M8 »Mg+l>Mg 


to 


ni+rt! 


"IMg +»lM g 


= 2.5, 


and for the II.C.2 “yin-yang I” 
orbit, the measured value = 3.20449 is closer to 

= 3.5. The only /3-orbit 


n+n 


nMg+«Mg 


= 3 than to 


n i+ni 

niM 8 +SlM 8 


whose measured value ij+- = 3.98059 is not changed by 
a change of word-reading puncture is the I.B.7 “drag¬ 
onfly”, as in both the “natural” and “permuted” cases 
= ni + Rl = 4 holds. Note also that for the 

nM 8 +^M 8 niMg+niMg 

I.B.6 “butterfly IV” orbit tjA- = 26.4324, which set the 
largest deviation (7.8%) with = 49/2 = 24.5, 


this deviation becomes twice as large with 
123/4 = 30.75. 


n i +m 

"lM g +"lMg 


T 


T s 



Figure 3: (color on line) The rescaled periods T r (w) of 
presently known (zero-angular-momentum) three-body orbits 
divided by the period of the figure-8 orbit T r (w s) versus one 
half the word length Ni w /2 in terms of “permuted puncture 
word reading”. Four (linear, quadratic, cubic and quartic) 
fits are shown as solid lines of different colors. 


2. Three digit (syzygies) sequences 

There is an alternative method of assigning a sequence 
of three symbols to any given “word” in F^. The rule 
for converting “words” written in terms of letters a, b, 


A, B to “numbers” in terms of three digits - (1, 2, 3) - 
is simple: Make the substitution a = 12, A = 21, b = 
32, B = 23 together with “cancelling in pairs rule”: 11 
= 22 = 33 = (empty sequence). So, for example abAB 
= (12)(32)(21)(23) = 12322123 = 123123 for the figure- 
eight, i.e. a an increase from N w = 4 to N s = 6 symbols. 
The number of symbols IV, changes to those shown in 
Table EO 


Table IV: Rescaled ratios of periods T of three-body orbits 
and Moore’s figure-8 period Tms, and ratios with period T^ 
of the first orbit ft in the given section of the Table, as func¬ 
tions of the numbers n a , n b , tia, n B , of a’s, b’s, A = a _1 ’s and 
B = b _1 ’s respectively, in the free-group word description of 
the orbit, as well as of the total length N a of the sequence of in¬ 
tegers 1,2,3 used to describe the topology, see text. Subscript 
P denotes the first, shortest orbits in a particular sequence of 
orbits. 


Label 

<T>2 

<t> 2 

n+n 

/ — \ 

N s 

N 

(Tm8>2 

< T S>2 

ng + fig 

(n, n) 

Ns(P) 

■L'S 

M8 

1 

1 

1 

it 

1 

6 

58 

0.999951 

0.999951 

1 

1,1 

1 

6 

I.B.l moth I 

2.62031 

1 

1 

2,3 

1 

18 

II.B.l yarn 

7.86391 

3.00114 

3 

6,9 

3 

54 

I.A.l butterfly I 

2.15774 

1 

1 

2,2 

1 

14 

I.A.2 butterfly II 

2.15762 

0.999944 

1 

2,2 

1 

14 

I.B.5 goggles 

4.29152 

1.9889 

2 

4,4 

2 

28 

I.B.7 dragonfly 

3.98059 

1 

1 

4,4 

1 

24 

I.A.3 bumblebee 

10.9534 

2.7517 

11/4 

11,11 

11/4 

66 

II.C.2 yin-yang I 

3.20449 

1 

1 

3,3 

1 

20 

II.C.3 yin-yang II 

12.8167 

3.9996 

4 

12,12 44/10 

88 

I.B.l moth I 

2.62031 

1 

1 

2,3 

1 

18 

I.B.3 butterfly III 

3.76742 

1.43778 

7/5 

3,4 

4/3 

24 

I.B.2 moth II 

4.63126 

1.76745 

9/5 

4,5 

14/9 

28 

I.B.4 moth III 

5.83013 

2.22498 

11/5 

5,6 

22/9 

44 

I.B.6 butterfly IV 

26.4324 

10.0875 

49/5 

24,25 

89/9 

178 


Note that many of the ratios N N ^ and in Ta¬ 

ble m are identical, but as many are not: in particular 
only trival ratios in the “moth I - butterfly IV” sub¬ 
sequence are identical. Perhaps more importantly, the 
ratios of N , which determine the ratios of periods 

, i.e., the slopes of linear sub-sequences are often 
not-as-good as with —+ 2 —: for example for “moth 
I” orbit, the measured value = 2.62031 is closer to 
—+?— = 2.5 than to , r +. , = 3. We leave it to the 

riM 8 +n.M 8 JV s (M 8 ; 

interested reader to compare the quality of the linear fit 
in Fig. 0] with that in Fig. |T] to N s and to draw their 
own conclusions. 
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T 


T s 



Figure 4: (color on line) The rescaled periods T r (w) of 
presently known (zero-angular-momentum) three-body or¬ 
bits divided by the period of the figure-8 orbit T r (wg) ver¬ 
sus the length N s of integer sequence divided by the length 
Ns = 6 of the figure-eight orbit integer sequence. Four (lin¬ 
ear, quadratic, cubic and quartic) fits are shown as solid lines 
of different colors. 
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